Molecular dynamics simulations of ethanol binding to the transmembrane domain of the glycine receptor: implications for the channel potentiation mechanism

The glycine receptor (GlyR) is potentiated by ethanol and other anesthetics. The potentiation mechanism at the molecular level is unknown and remains elusive, but mutagenic studies have shown that ethanol and other volatile anesthetics bind to a pocket between TM1, TM2, and TM3. The present study extends previous studies (Cheng et al., Proteins 2007;68:581-593) wherein we conducted homology modeling and molecular dynamics (MD) simulations to construct models of the homopentameric alpha1 subunits of the GlyR transmembrane domain in open and closed states. To understand the potentiation of GlyR by ethanol we compare the binding of ethanol molecules to the channel in these different states. We observe that ethanol stably resides inside solvent-accessible cavities found in the open state of GlyR that are formed by I229 (of TM1) in one subunit and S267 and A288 (of TM2 and TM3, respectively) in the adjacent subunit. The volume of these putative binding pockets is state-dependent. Selective binding to the open states of receptors has been proposed to explain the potentiating actions of this class of anesthetics. In accordance with this model, our MD simulations suggest that the potentiation of ethanol on GlyR may be effected through preferential binding of ethanol molecules to an inter-subunit binding pocket in the open state.     

Characterization of the functional role of the N-glycans in the AMPA receptor ligand-binding domain

The ligand-binding domains of AMPA receptor subunits carry two conserved N-glycosylation sites. In order to gain insight into the functional role of the corresponding N-glycans, we examined how the elimination of glycosylation at these sites (N407 and N414) affects the ligand-binding characteristics, structural stability, cell-surface expression, and channel properties of homomeric GluR-D (GluR4) receptor and its soluble ligand-binding domain (S1S2). GluR-D S1S2 protein expressed as a secreted protein in insect cells was found to be glycosylated at N407 and N414. No major differences in the ligand-binding properties were observed between the 'wild-type' S1S2 and non-glycosylated N407D/N414Q double mutant, or between S1S2 proteins expressed in the presence or absence of tunicamycin, an inhibitor of N-glycosylation. Purified glycosylated and non-glycosylated S1S2 proteins also showed similar thermostabilities as determined by CD spectroscopy. Full-length homomeric GluR-D receptor with N407D/N414Q mutation was expressed on the surface of HEK293 cells like the wild-type GluR-D. In outside-out patches, GluR-D and the N407D/N414Q mutant produced similar rapidly desensitizing current responses to glutamate and AMPA. We therefore report that the two conserved ligand-binding domain glycans do not play any major role in receptor-ligand interactions, do not impart a stabilizing effect on the ligand-binding domain, and are not critical for the formation and surface localization of homomeric GluR-D AMPA receptors in HEK293 cells.     

Predicted structure of the extracellular region of ligand-gated ion-channel receptors shows SH2-like and SH3-like domains forming the ligand-binding site.

Fast synaptic neurotransmission is mediated by ligand-gated ion-channel (LGIC) receptors, which include receptors for acetylcholine, serotonin, GABA, glycine, and glutamate. LGICs are pentamers with extracellular ligand-binding domains and form integral membrane ion channels that are selective for cations (acetylcholine and serotonin 5HT3 receptors) or anions (GABAA and glycine receptors and the invertebrate glutamate-binding chloride channel). They form a protein superfamily with no sequence similarity to any protein of known structure. Using a 1D-3D structure mapping approach, we have modeled the extracellular ligand-binding domain based on a significant match with the SH2 and SH3 domains of the biotin repressor structure. Refinement of the model based on knowledge of the large family of SH2 and SH3 structures, sequence alignments, and use of structure templates for loop building, allows the prediction of both monomer and pentamer models. These are consistent with medium-resolution electron microscopy structures and with experimental structure/function data from ligand-binding, antibody-binding, mutagenesis, protein-labeling and subunit-linking studies, and glycosylation sites. Also, the predicted polarity of the channel pore calculated from electrostatic potential maps of pentamer models of superfamily members is consistent with known ion selectivities. Using the glycine receptor alpha 1 subunit, which forms homopentamers, the monomeric and pentameric models define the agonist and antagonist (strychnine) binding sites to a deep crevice formed by an extended loop, which includes the invariant disulfide bridge, between the SH2 and SH3 domains. A detailed binding site for strychnine is reported that is in strong agreement with known structure/function data. A site for interaction of the extracellular ligand-binding domain with the activation of the M2 transmembrane helix is also suggested.     

Biotechnology of β-adrenergic receptors

The emergence of Biotechnology has provided pharmacologists with a variety of methods for investigating the structure, the function, and the regulation of membrane-bound receptors with a precision that was not imagined even five years ago. These new tools have been developed and used to analyze the known catecholamine β₁- and β₂ receptors and to discover and study a new subtype, the β-adrenergic receptor. We review here the salient features of each of these three receptors, compare their structural and functional properties, and propose models to explain their differential regulation in time and space. A whole family of proteins has now been found to share with the β-adrenergic receptors their most prominent features, including seven transmembrane domains and coupling with GTP-binding “G” proteins. We therefore propose that the biotechnology-based procedures developed for the β-adrenergic receptors will be well applicable to the other members of this “R₇G” family of receptors.     

Recombinant human PPAR-beta/delta ligand-binding domain is locked in an activated conformation by endogenous fatty acids

High-resolution crystallographic structures of recombinant human peroxisome proliferator-activated receptor ligand-binding domain (isotype beta/delta) reveal a fatty acid in the binding site. Mass spectrometry confirmed the presence of C16:0, C16:1, C18:0 and C18:1 in a ratio of approximately 3:2:1:4 with 11, Z-octadecenoic acid (cis-vaccenic acid) identified as the predominant species. These are endogenous fatty acids acquired from the bacterial expression system, and serve to lock the ligand-binding domain into the activated conformation. A requirement for crystal growth, the additive n-heptyl-beta-d-glucopyranoside, binds near the activation function helix where recognition of co-activator proteins occurs. Our observations suggest potential physiological ligands for human PPAR-beta/delta and highlight that reported binding studies must be treated with caution unless endogenous fatty acids have been removed from the sample prior to analysis.     

Structural rearrangements in the thyroid hormone receptor hinge domain and their putative role in the receptor function

The thyroid hormone receptor (TR) D-domain links the ligand-binding domain (LBD, EF-domain) to the DNA-binding domain (DBD, C-domain), but its structure, and even its existence as a functional unit, are controversial. The D domain is poorly conserved throughout the nuclear receptor family and was originally proposed to comprise an unfolded hinge that facilitates rotation between the LBD and the DBD. Previous TR LBD structures, however, have indicated that the true unstructured region is three to six amino acid residues long and that the D-domain N terminus folds into a short amphipathic alpha-helix (H0) contiguous with the DBD and that the C terminus of the D-domain comprises H1 and H2 of the LBD. Here, we solve structures of TR-LBDs in different crystal forms and show that the N terminus of the TRalpha D-domain can adopt two structures; it can either fold into an amphipathic helix that resembles TRbeta H0 or form an unstructured loop. H0 formation requires contacts with the AF-2 coactivator-binding groove of the neighboring TR LBD, which binds H0 sequences that resemble coactivator LXXLL motifs. Structural analysis of a liganded TR LBD with small angle X-ray scattering (SAXS) suggests that AF-2/H0 interactions mediate dimerization of this protein in solution. We propose that the TR D-domain has the potential to form functionally important extensions of the DBD and LBD or unfold to permit TRs to adapt to different DNA response elements. We also show that mutations of the D domain LXXLL-like motif indeed selectively inhibit TR interactions with an inverted palindromic response element (F2) in vitro and TR activity at this response element in cell-based transfection experiments.     

A theoretical analysis of HLA-DRbeta1*0301-CLIP complex using the first three multipolar moments of the electrostatic field

Interactions between the HLA-DRbeta1*0301 molecule and several occupying peptides obtained from computational substitutions made to the CLIP peptide are studied. The exploration was carried out using a vector composed of the first three terms of the multipolar expansion of the electrostatic field, namely, charge (q), dipole (d) and quadrupole (C). Comparisons between pocket-peptide interactions established that the binding pockets for this HLA molecule are ordered in terms of their importance for binding peptides, as follows: P1 > P4 > P6 > P7 > P9. A set of electrostatically distinct amino acids that determine interaction stability and specificity were identified for each pocket. The beta74R residue was especially identified as being the key amino acid mediating the occupying peptide binding for pocket 4; this residue has been recently associated with Graves' disease.     

Anatomy of protein pockets and cavities: measurement of binding site geometry and implications for ligand design.

Identification and size characterization of surface pockets and occluded cavities are initial steps in protein structure-based ligand design. A new program, CAST, for automatically locating and measuring protein pockets and cavities, is based on precise computational geometry methods, including alpha shape and discrete flow theory. CAST identifies and measures pockets and pocket mouth openings, as well as cavities. The program specifies the atoms lining pockets, pocket openings, and buried cavities; the volume and area of pockets and cavities; and the area and circumference of mouth openings. CAST analysis of over 100 proteins has been carried out; proteins examined include a set of 51 monomeric enzyme-ligand structures, several elastase-inhibitor complexes, the FK506 binding protein, 30 HIV-1 protease-inhibitor complexes, and a number of small and large protein inhibitors. Medium-sized globular proteins typically have 10-20 pockets/cavities. Most often, binding sites are pockets with 1-2 mouth openings; much less frequently they are cavities. Ligand binding pockets vary widely in size, most within the range 10(2)-10(3)A3. Statistical analysis reveals that the number of pockets and cavities is correlated with protein size, but there is no correlation between the size of the protein and the size of binding sites. Most frequently, the largest pocket/cavity is the active site, but there are a number of instructive exceptions. Ligand volume and binding site volume are somewhat correlated when binding site volume is < or =700 A3, but the ligand seldom occupies the entire site. Auxiliary pockets near the active site have been suggested as additional binding surface for designed ligands (Mattos C et al., 1994, Nat Struct Biol 1:55-58). Analysis of elastase-inhibitor complexes suggests that CAST can identify ancillary pockets suitable for recruitment in ligand design strategies. Analysis of the FK506 binding protein, and of compounds developed in SAR by NMR (Shuker SB et al., 1996, Science 274:1531-1534), indicates that CAST pocket computation may provide a priori identification of target proteins for linked-fragment design. CAST analysis of 30 HIV-1 protease-inhibitor complexes shows that the flexible active site pocket can vary over a range of 853-1,566 A3, and that there are two pockets near or adjoining the active site that may be recruited for ligand design.     

A large deletion/insertion-induced frameshift mutation of the androgen receptor gene in a family with a familial complete androgen insensitivity syndrome

Androgen insensitivity syndrome (AIS) is an X-linked recessive genetic disorder with a normal 46, XY karyotype caused by abnormality of the androgen receptor (AR) gene. One Chinese family consisting of the proband and 5 other members with complete androgen insensitivity syndrome (CAIS) was investigated. Mutation analysis by DNA sequencing on all 8 exons and flanking intron regions of the AR gene revealed a unique large deletion/insertion mutation in the family. A 287 bp deletion and 77 bp insertion (c.933_1219delins77) mutation at codon 312 resulted in a frameshift which caused a premature stop (p.Phe312Aspfs*7) of polypeptide formation. The proband's mother and grandmother were heterozygous for the mutant allele. The proband's father, uncle and grandfather have the normal allele. From the pedigree constructed from mutational analysis of the family, it is revealed that the probably pathogenic mutation comes from the maternal side.